61372 PIDHM I~ eh>er 6-71 


es 


FOREIGN TECHNOLOGY DIVISION 


AD ¢O9U4I90 


FOUR-LEVEL METHOD FOR PREDICTING GEOPOTENTIAL IN 
THE TROPOSHPERE AND LOWER STRATOSPHERE 
TAKING RADIATION AND TURBULENT 
HEAT TRANSFER INTO CONSIDERATION 


by 


Ween Llane eT ODOV. and Si. 1. Rudenko 


SATE WATER SURVEY LIVISIUR 
LIBRARY COPY 


npR 0 3 1974 


Approved for public release; 


| distribution unlimited. 


551.63 Galin, M.B. 

G1372 FOUR-LEVEL METHOD 

07110407 FOR PREDICTING 
GEOPOTENTIAL IN THE 
TROPOSPHERE AND 
LOWER STRATOSPHERE 


DEMCO 


ET DEM Teciterecen ; 


EDITED MACHINE TRANSLATION 


FTD-MT-24-2156-71 


FOUR-LEVEL METHOD FOR PREDICTING GEOPOTENTIAL IN 
THE TROPOSHPERE AND LOWER STRATOSPHERE TAKING 
RADIATION AND TURBULENT HEAT TRANSFER INTO 
CONSIDERATION 


Seen iat, tn hs ee LODOY <alldr oo, 2) hUudenko 
Fnglish pages: I17 


Bource: Voprmosy Gidrodinamicheskoy Teorii Kiimata 
i Dolgosrochnogo Prognoza Pogody 
(Questions Concerning the Hydrodynamic 
Theory of Climate and Long-Term Weather 
Pewecascine je Now leo) 900.,. 0D ae l= sue 

mequescver se PLD/PTR 

This document is a SYSTRAN machine aided 


translation, post-—-edited for technical ow 
Sccunacy Uys Ale Raul Wa Reltr ar. on 
BN 
Sp aS 
Rs 
\s oo 
Ag yy? woth 
> 9> 
ye 


Approved for public release; 
Gistributtonsuntimi Ted. 


THIS TRANSLATION IS A RENDITION OF THE ORIGI- 
NAL FOREIGN TEXT WITHOUT ANY ANALYTICAL OR 
EDITORIAL COMMENT. STATEMENTS OR THEORIES PREPARED BY: 
ADVOCATED OR IMPLIED ARE THOSE OF THE SOURCE 


AND DO NOT NECESSARILY REFLECT THE POSITION TRANSLATION DIVISION 
OR OPINION OF THE FOREIGN TECHNOLOGY DI- FOREIGN TECHNOLOGY DIVISION 
VISION. WP-AFB, OHIO. 


FID-MT- 24-22156-71 Date28 July19 72 


Digitized by the Internet Archive 
in 2023 with funding from 
University of Illinois Urbana-Champaign Alternates 


htips://archive.org/details/fourlevelmethodi0Ogali _ 


U. S. BOARD ON GEOGRAPHIC NAMES TRANSLITERATION SYSTEM 


Block Italic Transliteration Bioeck “italic [rans )iteration 
A a A a A, a F oP Comey Pog? as 

B 6 B 6 B, b Cenc Ce Sie te 

B Bs Be Mioats iat Tom ee 

Lair foe er es y y Vier y Dee 

pat | qq @ Deed > o @ Peas 

Ee Ee Yeiliyes- Fhe" xX ex xX tx Kh, kh 

K «K dK OK TA pA EV eu LEaG Loe cS 

3 3 fe Ci YA A Ca yy Chis acCn 

H uw Hu ioe ks Ww To te on, sh 

nA a Aa Yip Wr teed Ww Sheh, shch 
K kK K xX Re ae be? Dy tb " 

Sheen JI a ire a bl Al ow tan 

M ™ M « M, m b »b b eb ' 

haan H Le ag 39 #3 Fe 3 E, e 

Oo Oo Can OO »w tO »w Yuy «yu 
Moon Ton P, p A A AR Ya, ya 


SF Vevini tially, after, vowelserandeafter ny ss © elsewhere: 
When written as & in Russian, transliterate’ as yé or 4. 
The use of diacritical ma@rks’ is preferred, but such marks 
may be omitted when expediency dictates. 


FTD-MT-24-2156-71 ai 


FOLLOWING ARE THE CORRESPONDING RUSSIAN AND ENGLISH 


DESIGNATIONS OF THE TRIGONOMETRIC FUNCTIONS 


FTD-MT-24-2156-71 


Russian English 
sin sin 
cos cos 

tg tan 
ctg cot 
sec sec 
cosec csc 

sh sinh 
ch cosh 
th tanh 
cth coth 
sch sech 
csch esch 
arc sin sin71 
arc cos cos7+ 
arc tg tan-1 
arc ctg cot 
arc sec sec™ 
arc cosec cac- 
arc sh sinh71 
arc ch cosh71 
arc th tanh71 
arc cth coth~1 
arc sch sech~ 
arc csch esch7l 
rot curl 
lg log 


vB k 


FOUR-LEVEL METHOD FOR PREDICTING GEOPOTENTIAL 
IN THE TROPOSHPERE AND LOWER STRATOSPHERE 
TAKING RADIATION ‘AND TURBULENT HEAT 

TRANSFER INTO CONSIDERATION 


Weep Vir We Ks Popov “and S&S.) .. thudenko 


The linear four-level scheme for pre- 
dicting geopotential at isobaric surfaces 
au sone 80, 150, 250," anes 9e0 mb for’ the 
earth's northern hemisphere is used. A 
computational method is. proposed and examples 
are given for predicting (geopotential) for 
upecoy irec days. 


In the present work, a four-level model of baroclinic atmosphere 
ie used for predicting geoporentianl sy To-caiculate the inflow of 
heat from radiation and turbulent'neat transfer we use the method 
set forth in [1]. As our initial equations, we take the equation 
of vortex transfer sand the equation iof the inflow. of heat “in quasi- 


solenoidial approximation [2]: 


Pye 1 ow 
a + asin (Y» AY + 20a? cos 6) = — Qwa?cos 4 , (1) 


Pog se W 
Sepak he eo (2) 


¢ 


oT’ 1 ' 
nyse OPT hee tae 


where J is the flow function, T' - is deviation of temperature 
frou standard distribution, p (2) and P (2) - are standard distri- 
butions of density and pressure, t= P(z)/P(0), v: - is vertical 


velocTcy 
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Y ts-expressed. thus, i+ 


id 
1 2 02T’ P Pe » 


Z 
+e re-tde+ vet-*7,(2) (3) 


co 


where el oepade is optical thickness; 3 edhe is the density of the 


z 


substance which absorbs the radiation; 5 = is the coefficient of 
absorption of long-wave radiation: p=4roT3; WV’ - is the coefficient 
of turbulent heat transfer in a vertical direction; o - is the 
Stephan-Boltzman constant; C, - is the heat capacity of the air 


under constant pressure; Z=(€)_); 


We linearize equations (1) and (2) with respect to the zone 


flow, for which we present the functions in the following form: 


y= — a’a()cos6-+ 4”, 
T' =T,+ M(§)sin?6+-7", 
F=F-+F, 


where wp", T" and F" - are nonstationary nonzonal parts. 


From the relationship between the zonal parts of the functions 


of flow and temperature, T= — 4.20 coso » comes the relation 
S 
between parameter M and the index circulation a, 
20 » da 
MM sete ee Oe 


After linearization we obtain: 


) (] ” Oy” . 
(ar +2 Gr) 49 + 2(a +o) at = — 2a? cos 6 (4) 
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a eer 5 _ 
(ota 2 \ 77 — 2 cos Se ae ee — Man F", (5) 


As a boundary condition vertically we take 
W (0) = W (1)=0. 


We proceed in equations (4) and (5) to the finite differences 
Vertically. For this, we analyze the, @¢ntire thickness of the 


atmosphere with respect to variable § in 20 equal parts (such 


Brower = 0,05) by podntere; Eo: eye 519° Further, we write the 
equation of vortex transfer (4): for the four levels: E> E3> Es? 
19? where oy = 000.5: 53 =e Oy. b5 = 0.25 and E19 = 0.95. The 


equation for the inflow of heat (5) is written for three levels: 
Eos Eu. bio where E5 Te Vcges tate = Use. Ei =a eee UDC Le UG ane 
Porscvne derivatives ratios of Tinite differences, after several 
Cransiormations we obtain four equations fer wortex transfer and 
bares equations for the inflow of.iheat. Further, we use the 
condition of quasi-geostrophism and the barometric formula for 
replacing T". The derived equations are transferred according to 


(ae heinety we whl ret: 


% oy 
(ata YAR +4 (4 Seek Hl Past ey or pee 
—B, se (ME YS) Ha, (45 — 91) — 41,5 (95 — 95) = 0% (6) 


ies cif ih saucy acess mes) | ais 


7 
04, | ele: 
“Fs, rie Bs, Bane at %5) — By. 1 a (4; prs) ee a4 hes 1) Te 
+ a3 °5( 5 — 93) = 9; we. 
(Ge $55 Gr) [45 — TY = ¥5) + 1,25 (Hig — 94)] + 
0 Sy 
+ 85 9 id = P53 pT i (93+ $5) — Bs, — i auibsiee — 49) — 
— a, 5(%5;—4 3) = (8) 
Obi 
(Fr Tr %9 “5h 5) [Ate = a (%i9 —i?s)] + Bis, oS 
= Pio.s “ain e (95+ P19) =9; (9) 
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where 


aaa eee (45 be 
Ys 4w2a2g \ cos? 8 Jey 
- 1 ag—% , 
Bo = 2 (4, + %) — FF (%; — 4); Bs — Ba meek raat, ) ei 
Le 74. 
Give "ST Nree) Ase Wile 


B59 = 2(a;-+ 0) + (a, — 4) — (a; — 9)3 


___ 2 (a5 -— 43) |", 
Ps se er 


41,725 wt ie de 
3.1 Nr t 93.8 NE 


Bs,9 = 2 (25+) + (Gs — G3) = a (219 — 25); 


Gags . 0,62 (419 —a5) | 
85,3 = i ’ Bs, 19 = — Pig, sp 
3,9 _ 5,4 
oY eect Ya ye 45,3 > A 


Big, 9 = 2 (ag + 0) +> (a9 — 25), 


We search for the solution of the system (6)-(9) in the form 


de=C,e7 OF!" p™ (cogs), k=a1, 3, 5, 19. (10) 


To define the constant Cy we get a system of four algebraic 
equations. This system has a nontrivial solution if the determinant 
for the unknowns equals zero, from which the characteristic 
equation for the determination of o” 4s obtained. Cy and o™ are, 


generally speaking, complex, that 1S Cy=A,+iBr, o“=0,+i0;. 


The system indicated above of four algebraic equations permits 
the expression of the constants Cy for various levels by means of 
the analogous constants of any one level. In this instance 
we accepted the level corresponding to k = 5, which leads to the 


following relationships: 
A,= kh As crys RipBs, 
By =k .Bs + RinAs, 


pat BiG (11) 


where k and Key are constants, defining O- and O,: Since for 


k 
each n and m we have four values for o, then in place of relation- 
ships (11) we in fact have four relationships. Finally the 
solution takes the form: 


18 36 


b= i Y (Mi, cos mi. Nan Sin mh) P7 (cos 9), Grea) 


m=ln=m 


where 
Min = (At cos mart + Bz sin mort) fae + (A cos most + 
+. By sin mort) emit + (Az cos mort + By sin mart) emit ze 
++ (Ai cos mort + By sin m3;t) Bird \" : 
Nin = (Bt cos matt — Aksin malt) emt iy (Bz cos ma;t — 
— Ax sia ms;t) em + (Bz cos ms; — Aj sin ma; t) ese + 


4, ]m 
4 4 oe. 4 m3, t 
+ (Bzecos matt — Aj sin mo-t) e”™ ie 


feo OS arm ag he igo.) 


To determine AL and By we, use. Che initial.conditions: 


18 36 


yi — SS , ay cos mh + Be sin mh) P™ (cos 9), 


Mm =l fh= 1 


k=1, 3, 5, 19. (14) 


The sums on the right sides of (12) and’°(14) contain the 


terms corresponding to the even differences n - m. 


If in*relationships (13) we assume t = O' and use (11) and 
(14), we obtain for a determination of Ag and B, a system of 
equations of the eighth order, with coefficients which will be 
numbers k_, and k,,, but which in the right sides will have a? and 
BY Roi = ee nee es Oe 


We accomplish the transfer to geopotential using the relation- 


ship =n, where a - is the mean value of the Coriolis 
parameter. | | 


Computation is carried out according to the following scheme. 
We determine the initial expansion factors AY and BY, the mean- 


zonal values and circulation indexes a according 


17 ee bf -2 103 
to the initial geopotential fields. The standard functions - 
Legendre's polynomials Pe (cos 6) are considered each time and 
programmed in the machine's memory. Knowing the circulation 


indexes, we find the coefficients of the characteristic equation, 


0 0) 
K? By we find Ags Be 
by means of the solution’ of the eight-order system. Further, with 


the help of (11) we determine A, and B, for k = 1, 3, 19, while 


and then o. According to the values of o and A 


k k 
we find Men and Neen for any moment in time t according to formula 
(13). Summarizing series (12), we compute the desired values for 


the forecast fields. 


The initial fields of the altitudes of the isobaric surfaces 

150 and 250 mb were obtained by means of interpolation of the data 
of standard levels. Meanwhile it was assumed that the relationship 
between the altitudes of the isobaric surfaces were the same as 
for standard atmosphere. The altitudes of the isobaric surface 950 
mb were obtained at sea level pressure with the assistance of a 
barometric step, for which a value of 8.5 m/mb was accepted. In 
this way we used the following formulas: 

Hy 55 = 0,414F7 09 + 0,586 A 200, 

H1 955 = 0,456 A 199 77 0.5447 a99, 

F195) = 0,85 (py) — 950). 


Since Po is the sea level pressure in millibars, then for Hoso 


50? #150? Faso 
and Hoe is computed, we can recalculate the forecast values at 


we must get Values in decimeters. After forecast H 


standard levels according to formulas which are obtained similarly 
to oG1 5) 


HY 99 = 0,371 Ff sp + 0,629F7 25, 
Hag 22 048771, 25-7 0,563 a5), 
FA 399 = A125 — 120, 

Py = 950 + 1,176F7q;). 


The problem was programmed for calculation by an electronic 
computer. The forecast was computed for several days with the 
values for the geopotential of isobaric surfaces for each day 
printed. The calculation of the forecast for a 24 hour period 
for all four surfaces took-aboaut-10 minutes.. A-whole ‘series of 
examples was calculated, from which we reproduce two; one relates 
to the winter season (1 Feb 64), while the second relates to the 
summer (15 July 64). Zonal circulation in these two instances was 
substantially different. Thus, for 1 Feb. 1964 the zonal velocity 
monotonically increases with height, which is evident from the 
values of the circulation: indexes: M1 9/¥ Be ye Gul es os / uy = 0.067; 
O/W = 02072: a4 /w = 0.087. For 15“July 1964 the circulation has 
a Giiferent nature? Ateniecheairitodes..zbeginnine at the/50 mb 
level, we find an easterly zonal flow which is a general charac- 
teristic of the summer season. ~The clreulation indexes in this 
case take on the values: O14 g/W = 0.002; O15 /u = 0.047; O3/u = 0.039; 
a4 /w = 0.009. 


We examine the. forecast. actording to.the.initial data for 
1 Feb. 1964 (0300, Moscow time). A comparison of the forecast 
charts Ba50 fone iivie, 3-day perlods: with the wactual chants for 
that period shows that, as a whole, the synoptic situation was 
grasped by the forecast correctJly.. On the forecast chart for 2 
Feb. 1964 (the first 24 hours of the forecast) it is possible to 
novice a highespressure’ centers over. the Pyrenean Peninsula. On the 
actual map, 4. ridge corresponds ito this center. The position of 
troughs over Asia-Minor ,»the Himalayas, and the Pacific Océan were 
correctly given. In the second 24 hour period, high and low 
pressure centers over Central.Europe and the Mediterranean 
respectively were forecast accurately. In the third 24 hour period 
(Figures 1, 2) on the actual chart we see the weakening of the 
designated centers which is insufficiently reflected in the forecast. 


We examine the forecast charts AT for the first 24 hour period. 


200 
The position of troughs over the European territory of the USSR 
and the Atlantic is correctly. forecast. The ridge over the Urals 


was satisfactorily obtained, although with a slight displacement 
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0300 Moscow time, 4 Feb 1964, 


Chart ATS 0 a 
a) forecast (3 day); b) actual. 
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to the east. In Siberia, a low pressure in the region of 
Krasnoyarsk was correctly obtained. In the second 24 hour 

period, the troughs over the European territory of the USSR and 

a ridge over West Europe were also correctly forecast. The 
coincidence of forecast with the actual in the area of the Atlantic 
and the USA should be noted. In the third 24 hour period (Fig. 2) 
an area of low pressure over Siberia is displaced to the east on 
the forecast chart in comparison to its actual position. The 
situation over North America appears smoother than it actually 

was. The forecast chart of isobaric surface 950 mb was compared 
with the basis of the Surface pressure field. /in/the first 

24 hour period (Fig. 3), low pressure centers over Scandinavia, 
pressure centers over the eastern portion of the Mediterranean Sea, 
to the east of Japan, over the USA, including in the region of 

the Aleutian Islands, were forecast. In subsequent periods, 
centers over the westcoast of the USA, Japan, Scandinavia, and 
Siberia were satisfactorily obtained. The center over the Pyrenean 
was displaced to the east and exaggerated in comparison to the 
actual center. During the third 24 hour period, the forecast 

chart shows a high pressure region over the Pyrenean Peninsula, 

a low pressure region over the eastern portion of the Mediterranean 
Sea, and a vast high pressure area over Central Asia. All these 
areas coincide well with the analogous barometric formations on 

the actual chart. At the same time, two high pressure regions 
appear over the Pacific Ocean instead of the single one which 
appears on the actual map. The intense low in the region of the 


Labrador Peninsula appears weaker and is displaced to the west. 


We now examine an example of forecasting according to the 
source data for 15 July 1964. On the first 24 hour period, the 
formation of low pressure center at isobaric surface 50 mb is 
forecast on the ridge over the Atlantic, just as the formation 
of a high pressure center over the European territory of the USSR. 
The zonal circulation in the middle latitudes corresponds with 
the actual situation over the USA. In the second 24 hour period 


a high pressure center was discovered over Europe. The low 
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pressure centers over the Himalayas and Atlantic and the high 
pressure center over the Caspaian were exaggerated, In the third 
24 hour period, (Fig. 4) the low pressure center over the northern 
coast of Canada and a high pressure region over the USA and 
Scandinavia were forecast. On the 200 mb surface, high pressure 
centers were forecast over China, Arabia, Mexico; low pressure 
centers over the Okhotsk Sea, and troughs over the Atlantic and a 
ridge over the European territory of the USSR were also forecast. 
In the second 24 hour period the situation over Siberia, the far 
east of the USSR, and the USA was correctly depicted. The location 
of troughs over the Atlantic, the Black Sea, Central Asia, and 
Canada was satisfactorily predicted. In the third 24 hour period 
(Fig. 5) the situation in the middle latitudes was forecast better 
(the trough in the region of Greenland, the ridge over the northern 
Kuropean territory of the USSR, the trough over eastern Siberia 

and in the region of the west coast of the USA). The pressure 
centers in the southern latitudes appeared displaced to the east 
in comparison to their actual location. In the first 24 hour 
period on the 950 mb surface high pressure centers over Central 
Europe, the Atlantic Ocean, the western-most Pyrenean Peninsula, 
southern coast of the USA, north of Alaska and over the Pacific 
Ocean, and low pressure centers over Arabia, Tibet, eastern Siberia, 
and the west coast of the USA, (Fig. 6) were forecast. In the 
second 24 hour period, the meridional flow over Central Europe 

and the southwest flow over northwest Europe were forecast, as 

was the case in the first 24 hour period. The situation over the 
eastern USA and the.Pacific Ocean coincide with the actual. 
Cyclonic formations over Tibet, south of Greenland and over the 
Pacific Ocean east of Japan were reported. In the third 24 hour 
period (Fig. 7) a low pressure center over the Pacific Ocean was 
accurately forecast. At the same time, a low pressure belt, 

coming across the Asian continent, was split into two intensive 
centers on the forecast chart. The same applied to the high pressure 
region over the Atlantic. The high pressure region over Western 
Siberia appeared displaced to the northwest. 
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Fig. 4. Chart AT 
a) forecast (3 day); b) actual. 
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Fig. 7. Chart AT 
#8 July 1964, a)<forecast.<3day) 3b) 
actual. 
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